A Pre-Protostellar Core in L1551 
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ABSTRACT 

Large field surveys of NH3, C2S, "'^'^CO and C^^O in the L1551 dark cloud have revealed 
a prolate, pre-protostellar molecular core (L1551-MC) in a relatively quiescent región to the 
northwest of the weU-known IRS 5 source. The kinetic temperature is measured to be 9 K, the 
total mass is ^ ŽMq, and the average particle density is lO^'-lO^ cm~'^ . L1551-MC is 2Í25 x líll 
in projection oriented at a position angle of 133°. The turbulcnt motions are on the order of 
the sound speed in the médium and contain 4% of the gravitational energy, Egrav, of the core. 
The angular momentum vector is projected along the major axis of L1551-MC corresponding to 
a rotational energy of 2.5 x 10^'^ sin~^ i \ Egrav\- The thermal energy constitutes about a third of 
\Egrav \ and the virial mass is approximately equal to the total mass. L1551-MC is gravitationally 
bound and in the absence of strong, ~ 160 /íG, magnetic fields will Ukely contract on a ~ 0.3 Myr 
time scale. The hne profiles of many molecular species suggest that the cold quiescent interior 
is surrounded by a dynamic, perhaps infalling envelope which is embedded within the ambient 
molecular gas of L1551. 

Subject headings: ISM: clouds — ISM: individual(L1551) — stars: formation 



1. Introduction 

Stars play a vital role in the cosmos, enriching 
the interstellar médium and providing the basic 
building blocks for the formation of planets and 
life. Stars with masses < IMq are vastly more nu- 
merous than higher mass stars and contain most 
of the stellar mass in our Galaxy. Much advance- 
ment has been made toward the understanding of 
low mass star formation due in part to the study 
of nearby star-forming molecular complexes such 
as Taurus-Auriga. The dense molecular environ- 
ments most closely associated with the star forma- 
tion process have been identified (Beichman et al. 
1986; Myers 1987) and characterized in large sur- 
vey data (Jijina, Myers & Adams 1999, and ref- 
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erences therein). However, the initial conditions 
under which protostellar collapse occurs remain 
poorly understood. 

Dicarbon monosulphide is a powerful próbe of 
the dense interstellar médium in star-forming re- 
gions containing a multitude of emission lines in 
the millimeter and centimeter wavebands with no 
hyperfine structure (Saito et al. 1987), a large Zee- 
man splitting factor (Shinnaga & Yamamoto 2000) 
and high critical density (Wolkovitch et al. 1997). 
Observational and theoretical studies of G2S in the 
dense interstellar médium suggest that it traces 
the early evolutionary stages of pre-protostellar 
cores (Suzuki et al. 1992; Lee et al. 2003) and 
is otherwise not detectable in star forming re- 
gions (Lai & Crutcher 2000). 

In 2001, we conducted a C2S survey of Lynds 
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1551 (Lynds 1962; Snell 1981) in scarch of cnvi- 
ronments in the early stages of star formation. 
L1551 is an ~ SOM© dark cloud at a dištance 
of 140 pc (Kenyon, Dobrzycka & Hartmann 1994) 
associated with more than 20 known pre-main 
sequence stars (see Palla & Stahler 2002, and ref- 
ercnccs thcrcin). The encrgetics from thc cnibcd- 
ded young stellar objects are witnessed as molec- 
ular outflows (Snell, Loren & Plambeck 1980; 
Moriarty-Schicvcn, Butncr & Wannicr 1995), 
jets (Mundt et al. 1990) and numerous Herbig- 
Haro objects (Devine, Reipurth & Bally 1999). 

A single emission región appcarcxi in thc north- 
west quadrant of our C2S survey which we shall re- 
fer to as L1551-MC (Lynds 1551 Molecular Core). 
While L1551-MC is wcll separated from the known 
energetics in L1551, there is a low excitation 
Herbig-Haro object, HH265, which lies at the edge 
of the projected C2S emission (Devine, Reipurth 
& Bally 1999; Garnavich, Noriega-Crespo & Green 
1992). No known energetics in L1551 are directed 
toward this región (see above referenccs and also 
Figúre 1) and the driving source for HH265 re- 
mains a mystery. 

The many observations at centimeter and mil- 
limeter wavelengths uscd to investigate L1551-MC 
are described in §2 from which the morphologi- 
cal and physical properties are derived in §3. A 
discussion of thcsc rcsults is given in §4, and the 
conclusions, including future study plans, are sum- 
marized in §5. 

2. Observations and Reductions 

Table 1 displays the generál parameters of our 
observations while the text in this section sum- 
marizcs thc rclcvant observational Information not 
included in the Table. 

2.1. BIMA Observations 

2.1.1. X = 9mm 

On 2001 September 2 and 3 a - 25' x 25' re- 
gión centered on L1551 was mapped in C2S using 
the BIMA^ array in its most compact (D) con- 
figuration equipped with 30GHz receivers (Carl- 
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strom, Joy & Grcgo 1996). A 7 pointing mosaic 
(~ 10' x 10') centered on the peak emission de- 
tected in the D array coníiguration was then ob- 
scrvcd in the B and C array configurations in thc 
fall of 2001 and 2002. Only the C and D array data 
are used (see §3.3.2) which have sufficient u — v 
coverage to sample all spatial frequencies between 
27" and 4'.5. 

These observations utilized the highest resolu- 
tion corrclator configuration at BIMA (Wclch et 
al. 1996) which translates to a 0.054 km s~^ ve- 
locity resolution at the C2S transition frequency. 
Phasc calibration was cnsurcd with regular seans 
of the quasar 0449+113 while the flux scale of 
the observations was determined with seans of 
Jupiter and 3C84. All data wcrc rcduccd with 
the MIRIAD reduction package (Sault, Teuben & 
Wright 1995) using štandard MIRIAD tasks. 

2.1.2. X = 3mm 

The C2S emission región described above 
(L1551-MC) was then mapped at the transition 
frequencies of a numbcr of molecular traccrs in thc 
3 mm band. These mosaics all consist of a 7 point- 
ing, Nyquist sampled hexagonal grid centered on 
4h3r"09".l, +18°12'09'Í5 (J2000) covcring a - 2'.5 
región with 100% sensitivity. These tracks were 
run at the Hat Creek Rádio Observatory in the 
C and D array configurations bctwccn 2002 April 
and June. In addition to the spectral line Infor- 
mation, each correlator configuration contained 
a total of 800 MHz of continuum which was cali- 
brated and preserved. 

The u — v coverage was good at all molecular 
transition frequencies giving well sampled maps 
for sky angles 15" < 9 < 100" except for CS 
which was only observed in the C array configura- 
tion and is more sparsely sampled over this spatial 
range. The integration time per mosaic pointing 
totaled ~ 45 minutes for the CS map and ~ 2 
hours for all others. System temperatures were 
variable from ~ 250 to 650 K. All tracks contained 
5 minuté seans of 0530+135 or 0449+113 at in- 
tervals of < 25 minutes for phase calibration and 
5 minuté integrations on Saturn and W30II for 
flux calibration. All data were reduced with the 
MIRIAD reduction package. 



2 



2.2. GBT Observations of NH3 

Singlc pointing spcctra and spcctral linc maps 
containing the (1,1) and (2,2) inversion transition 
frequencies of the NH3 molecule were obtained 
with the NRAO^ 100 m Róbert C. Byrd Grecn 
Bank Telescope (GBT) in K-band on the nights 
of 2004 January 6th and lOth. The spectrometer 
was set to contain 50 MHz of bandwidth ccntered 
on a 23.7088144 GHz sky frequency using 9 lev- 
els of sampler resolution with 32768 lags in dual 
polarization mode. Thereforc, both the (1,1) and 
(2,2) inversion transition frequencies including the 
quadrupole sateUites fell within the spectral win- 
dow with 0.019 km spectral resolution. The 
absolute reference position was chosen to lie north 
of L1551 at 4'^31"20", +18°20'00" (J2000). 

The spectral line maps are the result of the 
On-the-Fly (OTF) mapping technique and cover 
approximately 5' x 5' centered on 4^31™08^.2, 
+18°13'07'Í5 (J2000). The GBT main beam ef- 
ficiency, r]B = 0.81 and the aperture eíHciency, 
riA = 0.58 in K-band. The 60arcsec s~^ scan 
rate and 13" row spacing used to make the maps 
sufBciently sample all spatial frequencies avail- 
able to the GBT within the mapped región. Sys- 
tem temperature calibrations performed every two 
rows minimize gain variations across the individ- 
ual maps and alternating the map scan directions 
sequentially between right ascension and decli- 
nation minimize systematic errors {e. g. striping). 
Typical systém temperatures were ~ 35 K on both 
obser ving nights. 

The single pointing spectra were calibrated and 
exported into ASCII tables with the AIPS++ Soft- 
ware package. All further reductions were per- 
formed using the IDL software of Research Sys- 
tems, Inc. The OTF maps were calibrated and 
gridded in AIPS-I-F then written to FITS formát 
for further processing in IDL. 

Our final spectra, including those in the spec- 
tral maps, are the combination of the two indepen- 
dent, calibrated polarization channels fit in off-line 
frequencies with a third-order polynomial. Our fi- 
nal spectral map is the combination by weighted 
mcan of six individual maps. Systematic effects 
from the OTF mapping technique are noticeable 

^The National Rádio Astronomy Observatory is a facility of 
the National Science Foundation, operated under coopera- 
tive agreement by Associated Universities, Inc. 



in the channels of the final map, although reduced 
significantly in the velocity integrated map. The 
total flux in the final spectral map was scaled up 
by 10% to match the total flux in our individual 
spectra taken at positions within the mapped re- 
gión. 

2.3. Kitt Peak 12 m Observations at A = 
3 mm 

2.3.1. On-the-Fly CO and C^^ O Maps 

A 20' x 20' región of L 1551 was mapped in 
1^00(1 0) on 2000 November 21 and 22 with 
the Kitt Peak 12 m Telescope^ using the OTF 
observing technique. The Millimeter Autocorre- 
lator (MAC) was configured in 2 IF mode with 
8192 channels and 150 MHz bandwidth. The main 
beam, taken to be a 2-dimensional Gaussian with a 
fuU width at half- maximum (FWHM) of 55", com- 
prises 90% of the fuU beam pattern. The aperture 
efSciency of the 12m at this observing frequency 
is 0.5. The 60arcscc s"-'^ scan rate and 19" row 
spacing of the maps provide spatial over-sampling. 
The number of rows per reference scan (OFF) and 
the number of OFFs per total power calibration 
varied over the course of the observations to keep 
the systém temperature variations below 3% be- 
tween measurements. System temperatures on the 
21st were very good, ~ 180 K, and decent on the 
22nd, ~ 275 K. 

A comparable OTF map in C^^0(1 0) ob- 
servcd on 2002 .lanuary 15 and on three consec- 
utive nights starting on 2002 October 21 covered 
the centrál 17' x 17' of L1551. The backend setup 
was identical to the ^^CO observations. System 
temperatures were good, ^ 250 K for our January 
observations and variable for the three nights in 
October. 

The raw OTF data were calibrated, converted 
to visibility data, baseline fit using a first order 
polynomial, gridded into an image and then writ- 
ten to FITS file formát with the AIPS reduction 
software. Further analyses were performed using 
IDL. 



'Tlic Kitt Peak 12 m Telescope is operated by tlie Arizona 
Rádio Observatory at Steward Observatory, University of 
Arizona, with partial support from the Research Corpora- 
tion. 
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2.3.2. CS and N2H+ Spectra 

Singlc pointing spectra wcrc takcn on thc night 
of 2002 November 30 at selected points within the 
C2S emission región at the transition frequencies 
of both N2H+(1 ^ 0) and CS(2 í). Thc fre- 
quency throw was set to ±1.0 MHz for CS and 
±3.0 MHz for N2H+. The MAC was configured 
in 2 IF mode with 12288 usable channels within a 
75 MHz band. 

Each spectrum was cahbrated and folded using 
UNIPOPS at the observatory and thcn output into 
FITS tables. A third order polynomial basehne fit, 
the combination of independent polarizations, and 
further analyses were all done using IDL. 

3. Results 
3.1. Overview 

Figúre 1 shows our fuU ^'^CO Kitt Peak 12 m 
map with known outflow sources overlaid in the 
Icft hand panel. Thc right hand panel is a closcup 
of the localized región of C2S emission L1551 re- 
vealed by our BIMA survey. The ^^CO is optically 
thick in this región and has a rclativc maximum 
which coincides with the peak of the C2S emission. 
HH265 lies at the eastern edge of the projected 
C2S. 

Figúre 2 shows a prolate, dense core seen in 
NH3 spatially ovcrlapping with the C2S emission. 
There are no compcUing young stcUar object can- 
didates near L1551-MC in either rádio, infrared, 
or X-ray data. Figúre 3 reveals the morphology of 
the NH3 emission to more closely rcscmble thc dis- 
tribution of C^®0 than thc optically thick ^^CO. 

Figúre 4 shows a comparison of many molecu- 
lar spectra taken at the peak of the C2S emission, 
dcfined as Point A, rclativc to thc centrál velocity 
of the NH3 emission. All lines in Figúre 4, exclud- 
ing C2S, are from the 12 m observations; BIMA 
observations of these same transitions did not re- 
sult in any detections, indicating the overall emis- 
sion is spatially smooth. The C2S emission peaks 
significantly to the red of the reference velocity. 
There are signs of self-absorption in both the ^^CO 
and CS spectra. Thc CS line profile in relation to 
the centrál velocity of the NH3 emission is con- 
sistent with an infalling envelope of L1551-MC. 
The N2H+ profile is slightly asymmetric, centered 
somewhat to the red of the NH3 centrál velocity, 



but is otherwise consistent with the NH3 emission 
profile. 

3.2. Is L1551-MC Starless? 

Jets, winds and outflows from young stellar 
sources can dramatically affect the interstellar en- 
vironment from which they are created. Therefore 
it is important to dctermine whethcr or not a stel- 
lar point source is present within a molecular core 
before the molecular line data are interpreted. 

The emission at all ffivlS' wavelengths in L1551 
is dominated by the L1551-IRS5/L1551NE and 
HL/XZ Tau regions and shows no evidence of any 
other sources within L1551. A recent paper by 
Gálfalk et al. (2004) reports the results from a 
deep mid-IR search for young stellar objects in 
L1551 using the Infrared Space Observatory (ISO) 
and four sources (sources 20, 23, 29 and 32 of their 
Table 2) lie within a 3' rádius of L1551-MC. Source 
32 is not considcrcd a potential protostar candi- 
date since it was only marginally dctected in the 
14.3 /xm band and has neither an optical nor near- 
IR counterpart. 

The thrce sources dctected at 6.7 //m are plot- 
ted on an (H — K) vs. {K — mg. 7) color-color di- 
agram (Figúre 4 in Gálfalk et al. 2004). Sources 
23 and 29 are consistent with dust reddened stars 
with no infrared cxcess. Source 20 was not de- 
tected in thc H band. However, this lower limit on 
niH places it above the reddening line at a position 
consistent with being a reddened main-sequence 
star. 

A SIMBAD search reveals one additional rádio 
source within the 3' rádius of the peak NH3 emis- 
sion, GRLll. This 1.4 GHz continuum source is 
most likely of extragalactic origin (Giovanardi et 
al. 2000; Sneh & Bally 1986). An X-ray source, 
BFR12, lies within 4" of GRLll and is likely re- 
lated (Bally, Feigelson & Reipurth 2003). An ad- 
ditional X-ray source from the Chandra survey of 
L1551, BFR2, is seen to the southwest of L1551- 
MC and is also most likely a background object. 

The combination of all our 3 mm BIMA con- 
tinuum data precludes the existence of a compact, 
< 10", source at ^ 100 GHz brightcr than 5.3 mJy 
with 99% confidence. Given these considerations 
we conclude that L1551-MC is starless. A plot of 
the Gálfalk et al. (2004), 2MASS and SIMBAD 
sources can be seen in Figúre 2. 
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3.3. Morphology 

3.3.1. Distribution and Orientation of Dense 
Gas 

The dense, inner regions of cores are traced 
well with NH3 emission due to a high abundance 
under these conditions and a tendency to avoid 
depletion onto grains (Aikawa, Ohashi &: Herbst 
2003). In this section the primáry morphologi- 
cal fcatures of L1551-MC arc characterized based 
on the NH3 emission in the región. Several sin- 
gle pointing spcctra taken at 23GHz within the 
larger L1551 cloud rcvcal that NH3 emission is 
not widespread but rather confined to the regions 
around L1551NE and IRS5 in addition to the re- 
gión mappcd in this papcr (scc Figúre 3). The 
results from these spectra are summarized in Ta- 
ble 2. 

L1551-MC appears roughly elliptical in projec- 
tion with an axial ratio, y = 2. Figúre 5 shows the 
velocity integrated NH3(1,1) emission in the opti- 
cally thin satellite componcnts as contours. This 
emission distribution is fit with a 2-dimensional 
Gaussian with the peak emission amphtude, 2-d 
centroid, major and minor widths, and position 
angle as free parameters. The gross morphologi- 
cal features, including the characteristic quantities 
from the average spectrum within a half-maximum 
NH3 contour are summarized in Table 3. 

A first moment map is created by fitting each 
spectrum in our NH3 data cube according to §A. 
Local velocity gradients calculated at each pixel 
by vector sum of the centrál velocity difference at 
the eight neighboring pixels in the íirst moment 
map arc shown in Figúre 5. There is some scatter 
to the local gradients, but there is an overall trend 
of higher velocities toward the southwest. Fitting 
a velocity gradient to the map in a manncr simi- 
lar to that done by Goodman et al. (1993) yields 

= 1.2 ± 0.1kms~^pc~^ at a position angle of 
224 ± 5°, aligned along the short axis of the veloc- 
ity integrated emission. The vector addition of all 
local velocity gradients in Figúre 5 yields an identi- 
cal gradient position angle but smaller amplitúde. 
These results suggest that L1551-MC is prolate. 
(We assume a prolate spheroidal morphology for 
L1551-MC in forthcoming calculations unless oth- 
erwise stated.) 

The velocity integrated emission in all satellite 
components averaged over elliptical annuli with an 



axial ratio of two ccntcrcd on the {a, S) coordi- 
nates in Table 3 at 13" rádius intervals is shown 
in Figúre 6. The 2-dimensional projected radial 
dištance from center is dcfincd as b, the impact 
parameter. The simplest characterization of the 
intensity profile is of a power law form, I{b) oc 
If this function is assumcd to dcscribc the intrin- 
sic shape of the intensity profile it becomes im- 
mediately apparent that it grossly overestimates 
the flux in the centrál pixels. A two component 
("broken") power law function of the form 



m 



C2b-P^ 



b<bo 
b>bo 



(1) 



where C2 = Ci6g^^ is much more succcssful 
in fitting these data. The fit according to the 
minimum chi-squared value of this intrinsic radial 
profile convolved with the GBT beam is shown in 
Figúre 6 as the solid line. There is little statistical 
uncertainty in our best fit parameters, but there 
is an unknown systematic error since the true 3- 
dimensional shape of L1551-MC cannot be deter- 
mined directly from our data. The štandard devia- 
tion in these parameters invoked by using elliptical 
annuli of varying axial ratios between y = 1 and 
y = 2 gives an estimate of this systematic error. 
The best fit parameter values arepi = 0.15 ±0.04, 
P2 = 1.48 ± 0.03, and bo = 0.040 ± 0.005 pc. 

An intrinsic intensity profile model of the form 



m 



1 + (V^o) ' 



(2) 



fits our data to even higher precision than the bro- 
ken power law. A similar model has bccn uscd 
by Tafalla et al. (2002) to fit the 3-dimensional 
radial profile of molecular cores with great suc- 
ccss and is a good approximation to a Bonncr- 
Ebert profile (Tafalla et al. 2004). This model 
is used here to fit the 2-dimensional intensity 
profile only for its convcnicncc and casily intcr- 
pretable parameters. The minimum chi-squared 
values for the free parameters of this function are 
lo = 2.4±0.2Kkms-i, feo = 0.05 ± 0.01 pc, and 
7 = 2.2 ± 0.1. This fit is displayed as the dashed 
line in Figúre 6. 

A broad brush picture of the mass distribution 
in L1551-MC can be inferred if it is assumed to 
be an optically thin prolate spheroid at constant 
temperature, constant abundance fraction, and in 
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local thcrmodynamic cquilibrium (LTE). The in- 
ferred fiat inner density profile with power index 
~ 1.2 and break at ~ 0.04 po where the power law 
index stecpcns to ^ 2.5 is typical for a starless 
core in Taurus {e. g. Caselli et al. 2002). 

This simple picture is most challengcd by the 
possibilities of a varying NH3 abundance fraction 
and a non-symmetric geometry in the third di- 
mension. Recent studies suggest that the abun- 
dance fraction of NH3 incrcascs toward the center 
of cores (Tafalla et al. 2004), so it is unlikely that 
an abundance fraction variation alone could pro- 
duce the flattened inner profile. It would also be 
unlikely that a geometrical effect could reproduce 
the efíect of having a centrally flat intensity dis- 
tribution. The existence of a turnover rádius at 
60 « 0.04 pc is therefore robust, while the slope 
of the profile both inside and outside bo is not. 
No significant N2H+ cmission is detected in the 
BIMA 3 mm data, sensitive to spatial scales be- 
tween 15 and 80" at a 3-0- level of O.lJybeam"^, 
suggesting that the densest gas in L1551-MC is 
distributed smoothly. 

3.3.2. The Envelope 

Carbon monoxidc, the sccond most abundant 
molecule in the Galactic interstellar médium, has 
a relatively small dipóle moment and traces low 
density molccular cnvironmcnts. Figurcs 1 and 3 
show ^'^CO and C^^O to be pervasive in L1551 and 
therefore these isotopomers próbe the low density 
envelope of L1551-MC as well as the environment 
in which it sits. 

If it is assumed that the rotational excitation 
temperature of ^^CO and C^^O are equal, the op- 
tical depths of these isotopomers can be estimated. 
The ratio of antenna temperatures of the ^^CO 
and C^^O transitions are 

ľI("CO) ^ (1-6-^3)^^3 

Tl(Ci80) (l-e---)r?i8 

where the subscripts 13 and 18 are used to de- 
note quantities related to ^^00 and C^*0 respec- 
tively and r] is the beam efRciency. The tele- 
scope parameters are nearly equivalent at these 
neighboring transition frequencies so to a good ap- 
proximation 7713 = jjig. It has been shown that 
for high column molecular gas {Ay > 3.5 mag) 
the ^^CO/C^*0 abundance ratio is relatively uni- 
form and consistent with terrestrial abundance 



ratios (McCutcheon et al. 1980) which implics 
ri3 = 5.4 Tig. Using these numbers, the mean op- 
tical depth of ^^CO in the NH3 emission región is 
1.9. 

The BIMA ^^CO and G^^O maps show a negli- 
gible amount of total emission. There is low level 
emission visible in the line wings of our ^"^CO map 
between the LSR velocities of 5.9-6.5 km s^^ and 
7.2-7.9 km s~^. However, there is not enough sig- 
nál abovc the 0.44 K per channel noisc; for any in- 
terpretation. No emission is detected in the line 
core where T13 > 1. The peak antenna tempera- 
ture in our Kitt Peak maps are 5.74 K and 2.23 K 
in "CO and C^^O respectively and our BIMA 
maps are well sampled between 15" and 80" with 
an RMS of ~ 0.4K. Therefore, the ^'''CO and 
C^^O are tracing a relatively smooth envelope of 
L1551-MC, with an insignificant amount of struc- 
ture smallcr than 80". 

The ^'^CO and C^^O emission extends smoothly 
from the región around L1551-MC into the greater 
L1551 cloud. The optically thin C^O follows 
more closely the overall form of NH3 emission 
showing no correlation to the C2S map. Therefore 
it appears that most of the mass in the envelope 
follows the gross morphology of the dense inner re- 
gions characterized in §3.3.1. This interpretation, 
however, does not account for depletion of these 
molecular species onto grains, which is known to 
occur in dense cores (Aikawa, Ohashi & Herbst 
2003; Lee et al. 2003; Tafalla et al. 2004). 

The C2S cmission in this región curls around 
the southwestern limb of the NH3 emission and 
stretches to the location of HH265 (see Figúre 2). 
The brightcst peak of C2S emission, dcfincd as 
Point A in this paper, coincides with a relative 
max;imum in the optically thick ^^CO emission 
(sce Figúre 1) but shows no such correlation in 
the optically thin C^^O. 

No C2S emission is seen to a 3-a level of 
0.27Jybeam~^ using BIMA in its B array con- 
figuration which samples spatial frequencies pri- 
marily between 10" and 35" at the our observing 
frequency. These data are not used in our final 
maps and spectra, yet the presented maps reflect 
an accurate morphology of this tracer given our 
scnsitivity limits. 

The CS spcctrum at Point A is liighly asym- 
metric containing a blue peak and a red slioulder. 
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consistent with self-absorption in an infalling laycr 
or shell (Myers et al. 1996; Evans 1999). The C2S 
emission is significantly redshifted with respect to 
thc systcmic NH3 vclocity. If thc C2S emission 
región lies on the near side of the core, the arcu- 
ate distribution of C2S may represent an infalhng 
laycr around L1551-MC. Thc possibihty of infall 
in L1551-MC is discussed further in §4.2. 

3.4. Physical Properties 

Many uscful quantitics dcscribing conditions 
within dense molecular environments can be de- 
rived from the NH3(1,1) and (2,2) inversion transi- 
tions (scc Ho & Towncs 1983). Our NH3(1,1) spcc- 
tra are fit according to the method described in 
§A. Minimum chi-squared values less than unity, 
as cxpcctcd for a dctailcd model, were typically 
achieved using this method and suggest that the 
LTE assumption is vaUd for this tracer. Table 4 
summarizes thc rcsults of our NH3 fitting and Fig- 
úre 7 shows fits to the spectra at Point A. 

3.4- 1- Kinetic Temperature 

Thc flux ratio of the NH3 (1,1) and (2,2) tran- 
sitions can be used to estimate the kinetic temper- 
ature of the emitting molecular gas (sce §A) . The 
bcst fit paramctcrs for thc (1,1) and (2,2) tran- 
sitions along with Equations A5 and A6 give a 
kinetic temperature at Point A, Tk = 8.7 ± 0.7 K. 

Thc scnsitivity of our NH3 OTF map was not 
sufficient to detect the (2,2) transition. While a 
variation in gas temperature by a factor of ~ 2 has 
been conjccturcd for Taurus corcs (Galii, Walms- 
ley & Gongalves 2002), the upper limit on from 
the Point B spectrum limits the temperature vari- 
ation to < 50% (sce Table 4). A constant kinetic 
temperature = 9 ± 1 K within L1551-MC is as- 
sumed in all following calculations which include 
temperature. 

3.4-2. Thermal and Turhulent Motions 

The full width at half-maximum thermal line 
width for a gas in LTE at kinetic temperature Tk 
is given by 

A«,\ = 81n2^ (4) 
m 

where m = 2.33 mjí is the mean particle mass of 
a gas of molecular hydrogen containing atomic hé- 
lium with a solar mass fraction {¥ = 0.28). Using 



Tk from §3.4.1, Avi/^ = 0.42 km s" ^ The non- 
thermal line width of a molecular species can be 
derived from the observed line width, Avobs, and 
kinetic temperature. 

Aí;^, = Aí;^,,-81n2-^ (5) 

whcrc rumoi is thc mass of the observed species. 
Assuming well mixed conditions, this quantity 
represents the non-thermal line width for the 
médium. The total line width can be obtained 
by adding the thermal and non-thermal compo- 
nents in quadrature. 

Thc mcan non-thermal line width is Au„t = 
0.20±0.002kms~^ and the mean total line width is 
Awťot = 0.46 ± 0.01 km s~-^ obtained by averaging 
the line widths from (^acli spectrum fit according to 
§A within a half-maximum NH3 contour. The er- 
rors in the line widths are from Landman, Roussel- 
Duprc & Tanigawa (1982) and by propagating the 
error in Tk. Thc sound speed in the médium 
is given by, (kTk/m)'^/'^ = 0.19 ± O.Olkms'S 
roughly cqual to Aw„í in thc región. Therc is no 
significant correlation between the projected diš- 
tance from the peak NH3 emission and the velocity 
width. 

3.4.3. Mass and Density 

The total mass of a molecular core cannot be 

observed dircctly and must be inferred using esti- 
mates of tracer abundances. The most abundant 
molecule, H2, does not radiate in the millimeter 
or centimeter bands duc to its homo-polar geome- 
try, leaving mass estimates of molecular regions 
dependent on conversion factors for less abun- 
dant moleculcs such as CO. Many processes af- 
fect molecular abundance ratios (see van Dishoeck 
& Blake 1998, for a comprehensive review) and 
published rcsults for X{CO) in Galactic molecu- 
lar clouds have varied by a factor of ~ 3 (Lacy 
et al. 1994). Following is an estimate of the total 
mass of L1551-MC first using our CO isotopomer 
data and NH3 data within a half-maximum con- 
tour of the velocity integrated NH3 intensity and 
then by consideration of the virial theorem. 

The total column of ^^CO in each pixel of 
L1551-MC is found assuming that ^^CO is ther- 
malizcd at 9 K. These values arc corrcctcd for op- 
tical depth using Equation 3 and then summed. 
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Using X(i3C0)= 1.12 x 10"^ converts this num- 
ber to a total mass of 1.9 M©. 

The total column of NH3 in the (1,1) state is 
derived approximating the partition function as 
the sum over all metastable states up to J = 
6 (Winnewisser, Churchwell & Walmsley 1979). 
The total NH3 mass of L1551-MC is found to 
be 2.4 x 10~^ Mq which implies a total molecu- 
lar mass between 0.8-2.4 M© for ^(NHs) between 
4.0 x 10"^ (Jijina, Myers & Adams 1999) and 
1.5 x 10-8 (Tafalla et al. 2004). The total mass of 
L1551-MC estimated from molecular line emission 
is Mmoi ^ 2 Mq good to within about a factor of 
2. 

The virial mass is obtained assuming virial 
equilibrium between thermal and gravitational 
self-energy. If it is assumed that L1551-MC is 
comprised of molecular hydrogen 



25 Aí;2 r 
24aGG ln2 



(6) 



Here ac is the geometrical and non-uniformity 

factor of Equation 7 which is of order unity (Bcrtoldi 
& McKee 1992), Av = Avtot from Table 5, 
and R is the semi-minor axis of the prolate 
spheroid. Equation 6 yields M^i^ = LQ^q'^Mq 
for aa = 0.87 

An estimate of the volume density from the 
NH3 fits is presented in Table 4. This density, 
riex , is estimated from the excitation temperature 
of the fit assuming non-LTE conditions and is a 
lowcr limit to the truc density of the médium 
in the derivation (see Equations A7 and foUow- 
ing). N2H+(1 - 0) is less likely to be in LTE than 
NH3(1,1) owing to a much highcr critical density. 
The hyperfine structure in the N2H"'"(1 — 0) tran- 
sition is fit in a similar fashion to the NH3 íits 
of §A only with the appropriatc centrál frcqucn- 
cies and line strengths for the 7 hyperfine compo- 
nents (Caselli, Myers & Thaddeus 1995; Tiné et 
al. 2000) input into Equation A4. The excitation 
temperature from the fit and temperature from 
§3.4.1 then determine the density of the médium 
using Equation A7 (scc also CascUi ct al. 2002). 
The best fit to the composite N2H+ spectrum gives 
an average particle density Ugx = 4 x 10* cm"^ in 
L1551-MC. 

The average density inside this core using a 
mass of 2Mq and volume of a prolate spheroid 
with dimensions from Table 3 is 5 x 10^ cm"^. In 



considcration of the factor of two crror in the mass, 
the unknown 3-dimensional geometry, and the es- 
timates of nex from the NH3 and N2H+ data, the 
average density in L1551-MC (n) « 10*-10^ cm"^. 

4. Discussion 

4.1. Energy Partition 

The gravitational stability and boundedness of 
L1551-MC can be estimated through investigation 
of the cnergy partition. 

The gravitational self-energy for an ellipsoidal 
mass with a power-law density distribution, p(r) oc 
r~P, is 

3GM2 



5 R 



where 



(1 — p I i) arcsine 
(1 - 2p/5) e 



(7) 



(8) 



and e is the eccentricity of the ellipsoid (Bertoldi 
& McKcc 1992). An estimate of the rotational 
energy can be derived by assuming a spherical ge- 
ometry with the rádius equal to the semi-major 
axis. 

2 

a,-Miž^„,^ (9) 



E, 



rot 



^mai . 2 • 
sm l 



where is the global vclocity gradicnt and aj = 
(1 — p/3)/(l — p/5). The thermal energy for a 
diatomic molecular gas is 5/2NkTk and the tur- 
bulcnt energy is estimated from the non-thermal 
component of the line width Eturb = 3/2M(7^(, 
where (T„í = 0.425Aí;„(. 

Table 6 summarizes the results for the energy 
partition in L1551-MC for threc diffcrcnt power 
law density profilcs, p = 0.0, 1.2, and 2.0, using 
both a spherical, y = í, and prolate spheroidal 
geometry with y = 2. In this table, Eyir = 2£^ + 
W and represents a measure of the gravitational 
stability based on virial argumcnts in the abscncc 
of magnetic fields (e.(?. McKcc & Zwcibel 1992). 
In this definition ^ = Etherm + Eturb + Erot is 
the kinetic energy term in the absence of pressure 
coníinement and W = Eg^av is the gravitational 
term. The negatíve values of E^i^ in Table 6 for all 
cases is strong evidence that L1551-MC is confined 
by its own self-gravity. The values most consistent 
with our analysis in §3 are p =1.2 and y = 2. 

The strength of the magnetic field is not known 
in this región. Given the total energy estimates 
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in Table 6, a field strength oí B > 160 /íG is 
needed for L1551-MC to be magnetically sub- 
critical which seems high based on Zeeman mea- 
surements (Lcvin ct al. 2001; Crutclicr, Hcilcs & 
Troland 2003), but is on the order of recent es- 
timates from polarized dust emission in starless 
corcs (Crutchcr ct al. 2004). The ambicnt niag- 
netic field position angle in L1551 was measured 
to be roughly 80° (Vrba, Strom, & Strom 1976) 
based on thc polarization of background stars in 
the periphery of the cloud. This angle is roughly 
consistent with the orientation of the jets and out- 
flows and may rcflect a relationship bctwecn the 
magnetic field geometry and the stellar formation 
process in L1551. The minor axis of L1551-MC 
is also in rough alignmcnt with thc ambicnt mag- 
netic field position angle, but further study of the 
magnetic field in the immediate vicinity of L1551- 
MC is nccdcd to dctermine conclusively its efíects 
on the evolution of this molecular core. 

To infer the force balance of a core from en- 
ergy considerations the cquilibrium state must be 
assumed. From §3.4.3, Mmoi ~ Myir suggest- 
ing that L1551-MC is virialized. While virial 
equilibrium is generally conjectured for ammo- 
nia cores (e.g. Jijina, Myers & Adams 1999), re- 
cent numerical results suggest that cores are not 
likely to be hydrostatic and that dynamic cores 
projected onto the sky can appear in equilib- 
rium (Ballesteros-Paredes, Klessen & Vázquez- 
Semadeni 2003, and references therein). Addi- 
tionally, there are no equilibrium conditions for an 
isolated, self-gravitating, prolate spheroid (Myers 
et al. 1991). Nevertheless, we can conclude from 
these analyses that L1551-MC is gravitationally 
bound and may be unstable. Thermal pressures 
providc some support against coUapse, but ulti- 
mately L1551-MC would have to be magnetically 
sub-critical to prevent collapse on a free-fall time 
scale, í// » 0.3 Myr. 

4.2. Core Environment 

L1551-MC appears fairly simple traced solely 
in NH3, but is significantly more complex when 
viewed in other molecular tracers. Figúre 4 shows 
the comparison between the molecular emission of 
sevcral tracers at Point A with thc vertical dot- 
ted line signifying the centrál LSR velocity of the 
NH3 (1,1) transition from Table 3. AU spectra in 
Figúre 4 are taken in the 3 mm band using the 



Kitt Peak 12 m Telescope except the C2S spec- 
trum which was taken from the BIMA 9 mm map, 
smoothed to ~ 55" to minimize beam dilution ef- 

fects. 

The ^'^CO spectrum at the top of Figúre 4 is op- 
tically thick; T13 = 1.9 ±0.2 using Equation 3 with 
the error determined by Monte Carlo simulation 
assuming no error in the isotopic ratio. The width 
of this line is much larger than expected from the 
direct conversion from the NH3 width and there is 
a dip in the line intensity at the centrál NH3 ve- 
locity. The ^^CO is likely tracing a column of gas 
which extends from the high deiisity interior into 
the more tenuous and turbulent surroundings. 

The CS line in Figúre 4 clearly shows a non- 
Gaussian shape which represents all of our CS 
spectra taken across the major axis of the C2S 
emission. This profile shape has a high peak at 
vlsr ~ 6.5kms~^ and a shoulder or second peak 
at vlsr ~ 7.0 km s"-'^ which could be due to an- 
other component along the line of sight or self- 
absorption. There is a ^ 50% decrement in in- 
tensity at the centrál velocity of the NH3 emis- 
sion, suggesting that this feature is due to self- 
absorption. The broad width of the CS line also 
means that this gas comprises a larger volume 
along the line of sight than the NH3. It is pos- 
sible that the asymmetric profile is due to gas in- 
fall from a layer (Myers et al. 1996) or spherical 
sheU (Evans 1999). 

Another intriguing feature seen in Figúre 4 is 
thc redshifted peak of thc C2S emission in rc- 
lation to the centrál velocity of the NH3. This 
0.16 km s~^ shift must be due to the kinematics 
of the gas in L1551-MC since the transition frc- 
quency used for our observations has an error (Ya- 
mamoto et al. 1990) much less than the observed 
line width. If the C2S emission lics on thc near 
side of L1551-MC, this implies that the C2S is 
moving toward the highest density environs and 
may bc further evidcncc of infall. Conversely, 
this kinematic shift could be due to outflow from 
L1551-MC. Chemical diíferentiation within molec- 
ular corcs is a well known phcnomenon (Tafalla et 
al. 2004, and references therein) and could also 
be contributing to the discrepancies seen between 
C2S and the other molecular tracers. Howcver, 
there is no clear evidence of depletion in either 
i^CO or C^^O. Alternatively, the C2S could be 
tracing a distinct component along the line of 
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sight. Givcn the suggestivc profile of the CS(2 — 1) 
line, the subsonic velocity shift of the C2S, and 
the coincidence of the C2S peak with the optically 
thick '^'^CO, the interpretation that the C2S gas is 
in a state of infall seems most likely. 

The N2H+ spectrum in Figúre 4 shows the iso- 
lated Jfi,f = lo,i ~^ (^1,2 component at Point 
A. The fit to the full spectrum (see §3.4.3) yields 
a value of tt = 5.5 ± 0.3 which means the op- 
tical depth of the isolated component is ~ 0.6. 
Assuming a beam filling factor of unity, the ex- 
citation temperature = 4.0 ± 0.1 K, and the 
total column of N2H+ within the 12 m beam is 
5 ± 0.5 x lO^^cm"^. The beam corrected col- 
umn depth ratio of our two nitrogenated molec- 
ular tracers, Af(NH3)/7V(N2H+)« 100, is consis- 
tent with other studies of starless cores in Tau- 
rus (Benson, CaselU & Myers 1998). The distri- 
bution of N2H+ emission peaks slightly to the red 
of the NH3 and is about a factor of 2 broader than 
that expected from 9 K gas. The N2H+ molecular 
emission is likely tracing a larger column of gas 
than the NH3. 

From these spectra, a picture emerges in which 
the cold, gravitationally bound inner región of 
L1551-MC is surrounded by a warmer, dynamic 
envelope showing evidence of gaseous infall. This 
envelope may not be uniform, but is connected to 
the more státie and smoothly distributed, tenuous 
molecular gas pervading L1551. 

5. Conclusions 

In this article, the diagnosis of a newly discov- 
ered molecular core found within the L1551 dark 
cloud is described. This core, referred to as L1551- 
MC, is located at 4'>3ľ°09''.9, +18°12'41" (J2000) 
and has a full width at half-maximum in NH3 
emission of ~ 2' x ľ at a position angle of 133° 
east of north. It has an inner región with a rel- 
atively fiat density profile which turns over to 
a more steeply declining function at a rádius of 
- 0.04 pc. L1551-MC has a mass of - 2Mq, 
an average density of lO^^-lO^ cm ^, and a kinetic 
temperature of 9 K which varies by less than 50% 
in the inner región. The velocity of the turbu- 
lence is on the order of the sound speed in the 
médium and comprises 4% of the gravitational en- 
ergy, Egrav, of the core. The rotational energy 
is ~ 2.5 x 10~^ sm~'^ i \Egrav\, while the thermal 



energy is 0.5\Egrav\- The core is gravitationally 
bound and either will coUapse on a time scale 
of 0.3 Myr or is supported by ~ 160 /xG mag- 
netic fields. L1551-MC has a high density inte- 
rior traced by NH3 and N2H+ suspended within a 
dynamic, perhaps infalling envelope which is con- 
nected to the more státie and tenuous molecular 
gas seen throughout L1551. 

From energy considerations it appears as 
though L1551-MC will contract to form a star or 
stellar systém. However, there is future work that 
could further characterize L1551-MC and more 
accurately plače it within the larger star forming 
context of L1551. An accurate measurement of 
the magnetic field in L1551-MC will further our 
understanding of the gravitational support mech- 
anism and constrain the core lifetime. The energy 
partition depends critically on the total mass and 
another independent method of determining the 
total mass, such as from dust emission, should be 
performed. While the state of the dense, inner 
regions of L1551-MC are presented herein, sup- 
plementary Information regarding the envelope of 
this core is needed. Spectral line maps of CS and 
its optically thin isotopomer, C"^^S, will determine 
the náture of the asymmetric CS line profile and 
provide a more complete physical, spatial, and 
chemical picture of the envelope. Lastly, HH265 
lies within the projected molecular emission of 
L1551-MC with no known driving source. Fur- 
ther investigation into the possible relationship 
between HH265 and L1551-MC may give deeper 
insight into the evolutionary state of L1551-MC. 

We thank Leo Blitz for allocating time for this 
project during the BIMA 1 cm observing season, 
and John Carlstrom for donating the usage of his 
30 GHz reccivcrs. J. J. S. would like to thank Carl 
Heiles for thoughtful input.This research has made 
use of the SIMBAD database, operated at CDS, 
Strasbourg, France. This research was partially 
supported by the National Science Foundation un- 
der grant AST 02-28963. Facilities: Berkeley- 
Illinois-Maryland-Association Array, Green Bank 
100 m Telescope, Kitt Peak 12 m Telescope 
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A. Ammonia Emission Model 



We model our observed spectra using 



TX{u)=A{l-e-^-^) 



where 
and 



hu 1 



MT) 



k exp(/ii//A;T) - 1 



For the (1,1) transition, the optical depth is 



18 

•E 



TT 2_^Sj exp 



{v 



2al 



(Al) 
(A2) 
(A3) 

(A4) 



where tt is the total optical depth of the transition, s j are the relative line strengths of the 18 hyperfine 
components, is the štandard deviation in frequency, and Vj are the centrál frequencies of the hyperfine 
transitions (Kukolich 1967; Rydbeck et al. 1977). To fit the (1,1) transition we wrote a wrapper for the 
IDL non-linear least-squares fitting routine MPFIT"' which easily allows for the tying together of parameters 
and returns all relevant fitting quantities such as the covariance matrix. More than 70% of the total (2,2) 
intensity lies in the main componont for which there arc 3 hyperfine transitions (Ai^i — O and AF = 0) 
with separations ~ IkHz (Kukolich 1967). Therefore, a single Gaussian component is used to fit our (2,2) 
spectra. 

Our non-linear least-squares fitting routine utilizes Equations Al and A4 to fit our (1,1) spectra with the 

amplitúde factor. A] the LSR velocity of the cmitting gas; the velocity width; and the total optical depth as 
free parameters. The velocity width is then analytically corrected for optical depth [e. g. Phillips et al. 1979) 
and instrumental broadening. The assumption that the excitation temperatures for all hyperfine transitions 
are equal is inherent in this fitting procedúre. 

The relative population of NH3 molecules in J = 2 and J = 1 states can be solved for assuming equal 
excitation temperatures and line widths for the (1,1) and (2,2) transitions (c.f. Ho et al. 1979) 



Tfl(2, 1) = -To 



In. 



-.911 Al 

922 A 



22 



V22 
Vil 



Tll S22 



In 



(A5) 



whcrc Tq is the energy difference between the J = 1 and J = 2 states in units of Kelvin, gn and (722 are the 
statistical weights, rn is the total optical depth of the (1,1) transition, and sn and S22 are the relative line 
strengths of the transitions. 

The kinctic tcmpcraturc, T^, can bo vcry wcU approximatcd from Tu if Tk ^ sincc transitions bctwccn 
K-ladders can only happen via coUision and non-metastable states within a K-ladder radiatively decay very 
quickly (Ho & Townes 1983; Walmsley & Ungerechts 1983; Stutzki & Winnewisser 1985). Then, only the 
(1,1), (2,2) and (2,1) states and thcir rcspective radiativc and collisional ratcs nccd to bc considcrcd in the 
detailed balance equation. Using Ho & Townes (1983) for the spontaneous emission rate and Danby et al. 
(1988) for the collisional rates, we obtain a transcendental equation for Tfc in terms of Tr and Tq 



Tfl = Tfc • 1 1 + ^ In [1 + 0.6 exp (-15.7/Tfc)]| 



(A6) 



'http; / /astrog.physics.wisc.edu/ ~craigm/idl/fitting.html 
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If a spccific valuc of thc bcam filling factor is assumcd, thc cxcitation tcmpcraturo can bc obtaincd from 
the amplitúde factor in our fits. The volume particle density can then be estimated using the equation of 
detailed balance for a 2-level systém (Stutzki & Winnewisser 1985; Caselli et al. 2002). 



T, 



íbg 



n 



crit {hv/k){l -T^lTk) 



where 



''crit 



P-—. 

Cul 



(A7) 



(A8) 



Aui and Cui are the spontaneous emission and cohisional de-excitation rates for the transition with u and / 
corresponding to the upper and lower energy levels within the (1,1) state. The escape probability, /? should 
be obtained by detailed modeling, but here we use the approximation, 



1-e- 



(A9) 



where r is taken to be the maximum optical depth from Equation A4. In Equation A7, 



1 — exp 



kT 



(AlO) 



Tk is the kinctic tcmpcraturc of thc mcdium, T^g = 2.73 K is thc background radiation temperature, and 
riex is the particle density of the coUisionally exciting species, taken to be H2. 
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4'^2™0' 45' 30° 15' 31™0' 45' 15' 12' 9' 6' 3' 31™0' 



Right Ascenslon (J2000) 

Fig. 1. — In the right panel, the velocity integrated ^^CO in L1551 is overlaid with pre-main sequence 
stars (stars), an incomplete sample of Herbig-Haro objects (triangles), known jet and outflow directions 
(arrows) (Mundt ct al. 1990; Devine, Reipurth & Bally 1999), and the spatial coverage of our NH3 map 
(dashed box). C2S emission is revealed to the northwest of center in a relatively quiescent región of L1551. 
In the right hand panel, a close-up view of the C2S emission región is shown overlaid on the ^^CO map which 
has been restrectched to the maximum and minimum values in the región. 
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Fig. 2. — Velocity integratcd NH3 cmission in greyscale with C2S contours at 30, 50, 70, and 90% the 
0.623Kknis~^ maximum ovcrlaid. Boxcs signify the point sourccs from Gálfalk et al. (2004), crosses are 
points sources from the 2MASS catalog, and the asterisks are SIMBAD objects. 
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4^31°'50^ 40^ 30^ 20^ 10^ 31™0' 

Right Ascension (J2000) 



Fig. 3. — Velocity integrated C^^O in L1551 with contours of NH3(1,1) overlaid at intervals of 30, 50 70 and 

90% thc pcak flux of 2.96 K km 8"^ The labclcd squarcs indicatc thc pointing ccntcrs of position switchcd 
23GHz spectra, the white star symbols are the pre-main sequence stars of Figúre 1, and the white triangle 
is the position of HH265. The Kitt Feak and GBT beams are shown from the bottom left respectively. 
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Fig. 4. — Spectra from an assortment of molecular transitions at Point A. The CS and N2H+ spectra are 
scaled by 5 x and the C2S spectrum is scaled by 2x. The vertical dotted hne signifies the systemic velocity 
of the NH3 emission at Point A. 
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Fig. 5. — First moment of NH3 emission in greyscale is overlaid with local velocity gradients. Arrow lengths 
signify the relative gradient strengths. The 224° position angle of the global 1.24kms^^pc^^ gradient is 
shown in the bottom right (length not to scalc). Velocity integrated NH3 intensity of the optically thin 
satellite components is shown in contour at 30, 50, 70 and 90% maximum intensity. 
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0.1 

b (pc) 

Fig. 6. — Intcgratcd intensity in the quadrupolc satcUitc componcnts of thc NH3(1,1) transition as a function 
of impact parameter. Error bars are the root mean square deviation of the mean intensity within an elhptical 
annulus and Ue mostly within the filled circles. Chi-squared is minimized for models derived from Equations 1 
and 2 and the results are plotted as the sohd and dashed lines respectively. 
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Fig. 7. — (a) NH3(1,1) spectrum at Point A with §A fit overlaid as the dashed line. The relative line 
strcngths of all 18 hypcrfinc componcnts arc dcpictcd as vcrtical lincs and thc LSR vclocity scalc in each 
quadrupole window is in reference to the strongest hyperfine line within. (b) The NH3(2,2) spectrum of 
Point A with single Gaussian fit shown as the dashed hne. Relative hyperfine line strengths are shown in 
the main component. 
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Table 1 

SUMMARY OF MOLECULAR LiNE DATA 



BIMA 



Transitions 


Frcqucncy (MHz) 


Au" (kHz) 


Obs. Type 


^heam 


C2S(32 ^ 2i) 


3375I.3737I 


6.1 


Mosaic 


27" x 20" 


i3C0(l ^0) 


110201.3541i 


24.4 


Mosaic 


10" x 9" 


Ci80(l ^0) 


109782.1734i 


24.4 


Mosaic 


10" x 9" 


CS(2 ^ 1) 


97980.9500^ 


12.2 


Mosaic 


9" x 7" 


N2H+(1 ^0) 


93176.2650^ 


12.2 


Mosaic 


11" x 10" 


C34S(2 ^ 1) 


109782.1734 


12.2 


Mosaic 


11" x 10" 


Green Bank 100m 


NH3(L1) 


23694.4955^'^ 


3.1 




34" 


NH3(1,1) 


23694.4955 


3.1 


APS'^ 


34" 


NH3(2,2) 


23722.6333^ 


3.1 


OTF 


34" 


NH3(2,2) 


23722.6333 


3.1 


APS 


34" 


Kitt Peak 12m 


i3C0(l ^ 0) 


110201.3541 


24.4 


OTF 


55" 


Ci«0(l ^0) 


109782.1734 


24.4 


OTF 


55" 


N2H+(1^0) 


93176.26502 


12.2 


FS^ 


65" 


CS(2 ^ 1) 


97980.9500 


12.2 


FS 


65" 



'^Post-reduction spectral resolution 
''On-the-Fly mapping 
■^Absolute Position Switching 
'^Frequency Switching 
iProm Pickett et al. (1998) 

^Prequency of the Jfi,f — lo,i ^1,2 hyperfine component (Caselli, 

Myers & Thaddeus 1995) 

^Hyperfine weighted mean transition frequency from Kukolich (1967) 

^Rydbeck et al. (1977) 
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Table 2 

SUMMARY OF THE GBT 23 GHZ POSITION SWITCHED ObSERVATIONS 



Source 


R.A.(J2000) 


Dec.(J2000) 


TB(l,ir 


Ts(2,2)^ 


<^Tb 


POINT A 


4:31:08.6 


+ 18:12:14.4 


1.41 


0.23 


0.04 


POINT B 


4:31:04.4 


+18:13:30.4 


1.09 




0.07 


IRS5 


4:31:34.1 


+18:08:04.1 


2.81 


0.53 


0.10 


L1551NE 


4:31:45.1 


+18:08:32.0 


0.44 




0.06 


XZTAU 


4:31:40.0 


+18:13:58.0 






0.06 


HH30 


4:31:37.5 


+18:12:26.0 






0.06 


OFFXZ 


4:31:35.8 


+18:14:45.0 






0.06 


MIDQNT 


4:31:25.0 


+18:13:30.0 






0.06 


BRTLMB 


4:31:15.5 


+18:06:55.0 






0.06 



Note. — Absence of data signifies a non-detection. 

'^Pcak value of main component brightness temperature in Kelvin. 
Beam fiUing factors are assumed to be unity. 
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Table 3 

Properties of NH:; Emission in L15r)l-MC 



R.A.'^ 


Dec.'' 


vlsr° 




maj.'^ 


min. 


P.A. 


(J2000) 


(J2000) 


(kms-1) 


(kms~^) 


(arcmin) 


(arcmin) 


degrees 


4:31:09.9 


+18:12:41 


6.72 


0.26 


2.25 


1.11 


133 



'^Coordinatcs of thc pcak emission. 

''The line profile characteristics are derived from the mean spectrum averaged 
within a half-maximum contour of the NH3 emission. 

'^"maj." and "min." are the FWHM of the major and minor axes respectively. 
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Table 4 

NH3 (1,1) AND (2,2) Fit Results 



Vf = 1 T^ = Tk 



Source 


(km S- M 


(kms-i) 


(K) 


TT 


T, 

(K) 


A^(NH.3) 

(cirr-) 


(cm-'') 


iV(NH3) 

(c-m-2) 


Vf 


POINT A 


6.72 


0.24 


8.7 


5.8 


4.4 


2.1 x 10^5 


5 x lO-"* 


3.7 x IQi'^ 


0.26 


POINT B 


6.71 


0.23 


< 12.4 


2.9 


4.5 


9.2 x lO^"* 


> 4 x 10^ 


> 1.5 x 10^'^ 


> 0.18 


IRS5 


6.42 


0.48 


10.0 


3.5 


6.4 


8.2 x lO^"* 


2 x 10-* 


4.1 x 10" 


0.47 


L1551NE 


6.55 


0.40 


< 20.6 


4.1 


3.2 


1.3 x 10^^ 


> 8 x 10^ 


> 4.6 x 10" 


> 0.02 



'^Uppcr limits to Tk dcrivcd from 3 <j Icvel in the (2,2) spectrum. Quantities dependent on were computed 
with uppcr limit valuc whcn appropriatc. 

^Density derived using equation Equations A7-A9 



Table 5 

L1551-MC Physical Characteristics'* 



Tk 


Mmol 




(n) 


^Vth 


^Vnt 


^Vtot 


(K) 






(cm-3) 


(kms-i) 


(kms-1) 


(kms-i) 


9 


2 


1.9 


5 x 10* 


0.42 


0.20 


0.46 



"^See text for a discussion of the uncertainties. 
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Table 6 

Energy Partition in L1551-MC 





y" 




OLI 


■^grav 


^therm 


Eturb X 10^ 


Erot^ XlO^ 


E <^ 










(10*3gi.gs) 


{\Egrav\) 


{\Egrav\) 


{\Egrav\) 


{\Egrav\) 


0.0 


1.0 


1.00 


1.00 


-0.9 


0.28 


3.2 


2.8 


-0.37 


1.2 


1.0 


1.15 


0.79 


-1.1 


0.24 


2.7 


1.9 


-0.46 


2.0 


1.0 


1.67 


0.56 


-1.5 


0.17 


1.9 


0.9 


-0.63 


0.0 


2.0 


0.76 


1.00 


-0.7 


0.37 


4.1 


3.7 


-0.18 


1.2 


2.0 


0.87 


0.79 


-0.8 


0.32 


3.6 


2.5 


-0.29 


2.0 


2.0 


1.27 


0.56 


-1.2 


0.22 


2.5 


1.2 


-0.51 



'^Power index for p{r) oc r~^. 

''Axial ratio, y = b/a 

^An inclination angle, i = 90° is assumed 
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